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We examine the general question of statistical changes experienced by ensembles of nonlinear 
random waves propagating in systems ruled by integrable equations. In our study that enters within 
the framework of integrable turbulence, we specifically focus on optical hber systems accurately 
described by the integrable one-dimensional nonlinear Schrodinger equation. We consider random 
complex fields having a gaussian statistics and an infinite extension at initial stage. We use numerical 
simulations with periodic boundary conditions and optical fiber experiments to investigate spectral 
and statistical changes experienced by nonlinear waves in focusing and in defocusing propagation 
regimes. As a result of nonlinear propagation, the power spectrum of the random wave broadens 
and takes exponential wings both in focusing and in defocusing regimes. Heavy-tailed deviations 
from gaussian statistics are observed in focusing regime while low-tailed deviations from gaussian 
statistics are observed in defocusing regime. After some transient evolution, the wave system is 
found to exhibit a statistically stationary state in which neither the probability density function 
of the wave field nor the spectrum change with the evolution variable. Separating fluctuations of 
small scale from fluctuations of large scale both in focusing and defocusing regime, we reveal the 
phenomenon of intermittency; i.e., small scales are characterized by large heavy-tailed deviations 
from Gaussian statistics, while the large ones are almost Gaussian. 


INTRODUCTION 


The field of modern nonlinear physics has sta rted with 
the pioneering work of Fermi and collaborators who 
studied a chain of coupled anharmonic oscillators, now 
known as the FPU system, with the aim of understand¬ 
ing the effect of nonlinearities in the process of ther- 
malization. Their unexpected results, i.e. the observa¬ 
tion of a recurrent behavior instead of the phenomenon 
of therm aliza tion, triggered the work by Zabusky and 
Kruskal 10l| who performed numerical simulations of 
the Korteweg de Vries equation (KdV), i.e. the long 
wave approximation of the FPU system, and made the 
fundamental discovery of solitons; such discovery lead to 
the development of a new field in mathematical physics 
that deals with integrable systems with an infinite num¬ 
ber of degrees of freedom. Som e years after the discov¬ 
ery, Zakahrov and Shabat |l02| found that the Nonlinear 
Schrodinger equation (ID-NLSE) is an integrable partial 
differential equations and has multi-solution just like the 
KdV equation. The following years were characterized 
by the search of new integrable equations and the study 
of their mathematical properties and solutions. 

In the late seventies and early eighties, besides soli¬ 
tons on a zero backgroun d, a new class of solutions of 
the ID-NLSE were found 103l4l06l| that describe the in¬ 
stability of a finite coherent background. Those solutions 
are s omet imes named “breath ers” : the classical Akhme- 
106l| and the Peregrine 1051 solutions breath just 


evolution variable and the perturbation of the coherent 
state is never small. 

The phenomenon of rogue waves (RW) in the ocean 
has been known to humans much before the discovery 
of integrability of partial differential equations; however, 
only in the last fifteen years a connection between those 
two fields has been made and it has been conjectured 
that the “breather” solutions of the focusing ID-NLSE 
could be considered as rogue wave prototypes |l07l llOSj . 
This idea has been rapidly picked up in diff erent scien¬ 
tific communities [5, ^,157, 8^ lid 111 | and a new 

field with fresh ideas and old equations has started. A 
first important step was the reproduction of the breather 
solutions of the ID-NLSE in water wave tanks [zllzl 
and in optical fibers [ 1 , 74, [t^. In order to generate 


diev 

once in their life and they describe in toto the modu- 
lational instability process ( also in its nonlinear stages). 
The Kuznetsov-Ma solution |l03l . Il04l| is periodic in the 


these coherent structures in controlled lab experiments, 
very specific and carefully-designed coherent initial con¬ 
ditions have been considered. However, in nature such 
conditions are almost never encountered; wind gene rates 
ocean waves via a non trivial mechanism ll.ll Ill4l | and 
the resulting wave field appears as a superposition of ran¬ 
dom waves characterized by Fourier spectra with a small, 
but finite, spectral bandwidth. This suggests that the 
problem of rogue waves should be investigated from a 
statistical point of view 33, M, Indeed, one of 

the major question to be answered in the field of rogue 
waves concerns the determination of the probability den¬ 
sity function (PDF) of the wave field for some given initial 
and boundary conditions. This is definitely not an easy 
task and, nowadays, given a nonlinear partial differential 
equation, there is no systematic theory that allows one 
to determine the PDF of the wave field. 
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The field of rogue waves has “belong” to oceanogra¬ 
phers until the the pioneering experiments with opti¬ 
cal fibers described in [23|. Since then, optical rogue 
waves have been studied in various contexts such as su¬ 
percontinuum generation in fibers d, [ 2 ^ 2^, 10, H, 7^ , 
propagation in optical fiber described by the “pure” ID - 
NLSE [iOl or with higher order dispersion 117, 118 |. 


M El 110, 111 and Raman fiber amplifiers |36l|. 


laser filamentation [10], passive cavities Isj. llTj . lasers 


From the general point of view and beyond the ques¬ 
tion of rogue waves, the field of nonlinear optics has 
then grown as a favorable laboratory to investigate both 
statistical properties of nonlinear random waves and 
hyd rodynamic-like phenomena [3, El m Ei, lilini 
[l20|- Indeed, the field of incoherent dispersiye wayes re¬ 
semble yery much the classical field of fluid turbulence 
where, instead of wayes, eddies interact with each other, 
giying birth to new eddies of different size. This mecha¬ 
nism is at the origin of the celebrated Kolmogoroy cas¬ 
cade of the three-dimensional turbulence which is char¬ 
acterized by a constant flux of energy within the so called 
inertial range. A source and a sink of energy are required 
in order to maintain the cascade. Many years after such 
concept was deyeloped, it was found that the cascade 
is intermittent, i.e. the statistical properties of the ye- 
locity field yary with the scales, becoming less Normal 
for smaller scales, see for example [sl for references. In 
the light of the paper Q, this idea will be discussed in 
the present paper in the context of the dynamics of in¬ 
coherent wayes ruled by the integrable ID-NLSE. Such 
equation proyides a bridg e between nonlinear optics and 
hydrodynamics, see |11()| for a one to one comparison. 
In particular, the focusing ID-NLSE describes at leading 
order the physics of deep-water waye trains and it plays a 


central role in the study of rogue wayes [i,|3l|,[5i,[6i[2i 
Moreoyer, the focusing ID-NLSE is the simplest partial 
differential equation that describes the modulational in¬ 
stability phenomenon that is belieyed to be a fundamen¬ 
tal mechanism for the formation of RW [s^, ■ 

As mentioned, such wayes emerge in the ocean from the 
interplay of incoherent wayes in turbulent systems. The 
theoretical framework combining a statistical approach 
of random wayes together with the property of integra- 
bility of the ID-NLSE is known as integrable turbulence. 
This emerging field of research recently introduced by 
V. Zakharoy relies on the analysis of complex phenom¬ 
ena found in nonlinear random wayes systems described 
by an integrable equation d, 0, QlilE^, EH- strictly 
speaking, the word “turbulence” is not fully appropriate 
in the sense that the dynamics in Fourier space is not 
characterized by a constant flux of a conseryed quantity 
because the system is Hamiltonian (no forcing and dissi¬ 
pation are included). For these integrable systems, giyen 
an initial condition, the spectrum generally relaxes to a 
statistically stationary state that in general is different 
from the standard thermal equilibrium characterized by 


the equipartition of energy. The prediction of the spectra 
of such final state and its statistical properties is the ob- 
jectiye of the integrable turbulence field. In the weakly 
nonlinear regime [42j | , starting with incoherent initial con¬ 
ditions in the ID-NLSE, deyiation from gaussian statis¬ 
tics has been predicted. In hydrodynamical numerical 
simulations performed with enyelope equations and ex¬ 
periments made in water tanks, non gaussian statistics 
of the waye height has also been foun d to emerge from 
random initial conditions 68, 6^, II5 |. 

While in the water waye context the NLSE is only 
a crude (but reasonable) approximation of the original 
equations of motion, the field of nonlinear fiber optics is 
a promising field for the inyestigation of integrable turbu¬ 
lence because optical tabletop “model experiments” ac¬ 
curately described by the ID-NLSE can be performed 
0, 0, El, Eli- Despite the numerous works deyoted 
to optical RW, the generation of extreme eyents from 
purely stochastic initial conditions in focusing ID-NLSE 
model experiments remains a crucial and open question 

BSEilel. 


In this paper, we reyiew and extend a number of re¬ 
sults recently obtained by the authors of this paper from 
optical fiber experiments 0 in the anomalous and 
normal dispersion regime. The dynamics of the wayes in 
the considered fiber is described with high accuracy by 
the focusing and defocusing ID-NLSE. In the focusing 
regime, the idea is to implement optical fiber experiments 
conceptually analogous to the water tank experiment de¬ 
scribed in [691 where wayes with a finite spectral band¬ 
width and random phases are generated at one end of the 
tank and the eyolution of the statistical properties of the 
waye field is followed along the flume. Using an original 
setup to oyercome bandwidth limitations of usual detec¬ 
tors, we eyidence strong distortions of the statistics of 
nonlinear random light characterizing the occurrence of 
optical rogue wayes in integrable turbulence. 

In the defocusing regime, modulational instability is 
not possible and the eyolution of incoherent wayes does 
not lead to the formation of rogue wayes. The statistics 
of waye intensity, initially following the central limit the¬ 
orem, changes along the fiber resulting in a decrease of 
the tails of the PDF. This implies that the probability 
of finding a rogue waye is lower than the one described 
by linear theory. Implementing an optical filtering tech¬ 
nique, we also report on the statistics of intensity of light 
fluctuations on different scales and we obserye that the 
PDF of the waye intensity show tails that strongly de¬ 
pend on the scales. This reyeals the phenomenon of 
intermittency, preyiously mentioned, that is similar to 
the one reported in seyeral other waye systems, though 
fundamentally far from being described by an integrable 
waye equation. We also report new and original experi¬ 
mental results for partially coherent wayes haying a broad 
spectrum. We demonstrate in particular the emergence 
of strongly non gaussian statistics with low tailed PDF 
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in the defocusing regime. 

The paper is organized as follows: numerical simula¬ 
tions of focusing and defocusing ID-NLSE are first con¬ 
sidered and described in Sec. □ Optical fiber experiments 
designed to investigate changes in the statistics of ran¬ 
dom light fields and showing results in agreement with 
simulations are presented in Sec. □ In Sec. „ we show that 
the integrable wave system under consideration exhibit a 
phenomenon of intermittency both in focusing and defo¬ 
cusing regime. In Sec. ^ we summarize our work and we 
discuss open questions about integrable turbulence. 


SPATIO-TEMPORAL, SPECTRAL AND 
STATISTICAL FEATURES ARISING FROM 
NONLINEAR PROPAGATION OF RANDOM 
WAVES IN SYSTEMS RULED BY THE 
INTEGRABLE ID-NLSE 


General framework and description of the random 
initial condition 


Our study enters within the general framework of the 
integrable ID-NLSE: 


ii>t + i^xx -f 2cr I'l/'PV' = 0 


( 1 ) 


where tjj{x,t) is the complex wave envelope. The pa¬ 
rameter a determines the focusing {a = -1-1) or de¬ 
focusing (cr = —1) nature of the propagation regime. 
In nonlinear fiber optics, it is relatively easy to ex¬ 
plore each of the two propagation regimes just by 
changing either the fiber or the wavelength of light 
[ 3 ^. Eq. (HJ) conserves the energy (or Hamiltonian) 
H = Hl + Hnl that has a nonlinear contribution 
Hjsil = —cr f \'ij}{x^t)\^dx and a linear (kinetic) contribu¬ 
tion Hl = f k'^\^{k,t)\‘^dk, the Eourier transform being 
defined as ili{k,t) = l/v^/ 'ilio{x,t)e~^^^dx. Eq. ([T]) 
also conserves the number of particules (or power) N = 
f \'ilj{x,t)\'^dx and the momentum P = J k\ip{k,t)\'^dk. 

In our paper, the initial conditions that are used are 
non-decaying random complex fields. We examine a situ¬ 
ation that is very different from the problem of the Fraun- 
hoffer diffraction of nonlinear spatially incoherent waves 
already considered in ref. 13 , 1 ^ 1^ • In these pa¬ 


pers, the nonlinear propagation of a speckle pattern of 
limited and finite spatial extension is studied in focus¬ 
ing and defocusing media. On the other hand we con¬ 
sider here continuous random waves of infinite spatial 
extension. This corresponds for instance to an experi¬ 
mental situation in which a partially coherent and con¬ 
tinuous (i.e. not pulsed) light source of high pow er is 
launched inside a single-mode optical fiber [3|, ISOf. In 
numerical simulations, the random waves are confined 
in a box of size L and periodic boundary conditions 
{ip{x = 0,t) =Mx = L,t)) are used to describe their 
time evolution ISOll. 


The random complex field 'il){x,t = 0) = ipoix) used 
as initial condition in this paper is made from a discrete 
sum of Fourier components : 

•0(0:,^ = 0) ='0o(a:) = ^ (2) 


with 'f/jQn = 1/T Jo '^oix)e ^^’^°^dx and feo = 27r/L. The 
Fourier modes ipon = are complex variables. 

In the random phase and amplitude (RPA) model, gen¬ 
eration of a random initial complex field is achieved by 
taking \ipon\ and (f>on as randomly-distributed variables 
( 5 ^ . Here, we will mainly use the so-called random 
phase (RP) model in which only the phases (pon of the 
Fourier modes are considered as being random [^. In 
this model, the phase of each Fourier mode is randomly 
and uniformly distributed between —tt and tt. Moreover, 
the phases of separate Fourier modes are not correlated 
so that < 6"^“™ >= 6nm- In the previous expres¬ 

sion, the brackets represent an average operation made 
over an ensemble of many realizations of the random pro¬ 
cess. 5nm is the Kronecker symbol defined by Snm = I if 
n = m and 6nm = 0 if n ^ to. 

With the assumptions of the RP model above de¬ 
scribed, the statistics of the initial field is homogeneous, 
which means that all statistical moments of the initial 
complex field ipoix) do not depend on x [13, Si- The 
power spectrum no(fc„) of the random field ipo{x) then 
reads as : 


^ '0On'0Om ^— riOn ^nm — rioikjP). 


(3) 


with kn = nkg. In the limit where L —>■ 00 , the fre¬ 
quency separation between two neighboring frequency 
components kn and kn+i tends to zero and the discrete 
spectrum no(fc„) becomes a continuous spectrum no(fc). 

The RP model is often used in the contexts of hydrody¬ 
namics where the power spectrum nojk) is given by the 
so-called JONSWAP spectrum [t^, Q also been 

used in optics where simple gaussian or sech profiles are 
often used for the function no{k) idlllS. 

In this Section, we consider a random complex initial 
field having a gaussian optical power spectrum that reads 


no(k) = no exp 






VAP ) 


(4) 


where A/c is the half width at 1/e of the power spec¬ 
trum. Fig. 1 shows a typical example of a partially co¬ 
herent complex field generated using the RP model and 
a power spectrum given by Eq. ®. Fig. 1(b) shows 
that the values of the spectral phases (po{k) are randomly 
distributed between —tt and tt. Fig. 1(c) shows the ran¬ 
dom evolution of the real part Rq{x) = ^{'tlJo{x)) of the 
initial field that is computed from the spectra shown in 
Fig. 1(a) and 1(b). We do not present here the evolu¬ 
tion of the imaginary part loix) = ^{tpo{x)) of the initial 
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Figure 1. (a) Gaussian power spectrum no(fc) of the initial 

condition such as defined by Eq. Q ( Ak — 0.5, no = 1.129). 
(b) Spectral distribution of the phases <j)oik) of the Fourier 
modes used to compute the random initial field ipo{x). (c) 
Real part Ro{x) = Sfi('ifio(x)) and (d) power Po(x) = \'4)o{x)\^ 
of the random field computed from spectra shown in (a), (b). 
(e) PDF of Ro{x) and (f) PDF of Po{x) showing that the 
initial condition has a gaussian statistics. 

field because it is qualitatively very comparable to what 
is shown in Fig. 1(c). However it is important to no¬ 
tice that the RP model produces a random field having 
real and imaginary parts that are statistically indepen¬ 
dent, i.e. < Ro{x)Io{x) >= 0 yX. The spatial evolution 
of the power Po{x) = \'ipo{x)\'^ = Rq{x) + I§{x) of the 
random initial field is shown in Fig. 1(d). Note that 
the numerical values of the parameters ng and Afe have 
been chosen in such a way that the number of particules 
N = 1/L |^(x,f)pdx is equal to unity. 

As previously described, the random complex field 
ipoix) used as initial condition is produced from the linear 
superposition of a large number of independent Fourier 
modes having randomly distributed phases. As stated 
by the central limit theorem, the statistics of the random 
process produced from such a superposition follows the 
normal law. The RP model thus produces a complex field 
having quadratures that are statistically independent and 
that have the same gaussian statistics. It is straightfor¬ 
ward to prove that the statistics of power fluctuations 
Pq follows the exponential distribution Sim. Note that 
the PDF for the fluctuations of the amplitude Aq{x) = 
\iPq(x)\ is given by the Rayleigh distribution defined by 
PDF[Aq/ < Aq >] = Aq/ < Aq > .exp{—AQj < Aq >) 

ilH 

In order to illustrate these statistical features from 
numerical simulations, we have performed the analysis 
of the statistical properties of the complex field gener¬ 
ated from the RP model by producing an ensemble of 


10“^ realizations of the random initial field. From this 
ensemble, it is in particular possible to compute the 
PDF for the fluctuations of Ro{x) = ^{tpo{x)) and of 
Poix) = |^/'o(a:)p. Fig. 1(e) shows that the PDF of 
Ro{x)/ < Ro{x) > is gaussian. The PDF of the normal¬ 
ized imaginary part Iq {x) / < Iq (x) > of the complex field 
is not shown here but it is rigorously identical to the PDF 
of Ro{x)/ < Ro{x) >. As it is illustrated in Fig. 1(f), 
the PDF of the power is given by the exponential distri¬ 
bution, i.e. PDF[Po/ < Pq >] = exp(—Pg/ < Po >)■ 

Focusing regime 

In this Section, we consider the focusing regime (cr = 
-|-I) and we use numerical simulations of Eq. ([1} to inves¬ 
tigate the propagation of a partially coherent wave gen¬ 
erated at t = 0 by the random complex field described 
in Sec. □ We will consider the spatio-temporal dynamics 
of the partially coherent wave and we will also discuss 
spectral and statistical changes occurring in time. 

Our numerical simulations have been performed by 
using a pseudo-spectral method working with a step- 
adaptative algorithm permitting to reach a specified level 
of numerical accuracy. The numerical simulations are 
performed by using a box of size L = 257.36 that has 
been discretized by using 4096 points. Statistical proper¬ 
ties of the random wave are computed from an ensemble 
of 10^ realizations of the random initial condition. 

Fig. 2(a) shows the spatio-temporal evolution of the 
partially-coherent wave seeded by a random initial field 
having properties that are described in SeCo and that are 
synthesized in Fig. I. At the initial stage of the nonlin¬ 
ear evolution (t ~ 0), the power |^(a::,t)p of the wave is 
slowly and randomly modulated. The spatial scale of the 
fluctuations of |^(x, t) p at t ~ 0 is determined by Ak and 
it is typically around 10 in the simulations shown in Fig. 
2 (see also Fig.1(d)). Numerical simulations show that 
series of peaks emerge from the random initial condition. 
The density of these peaks is higher is those regions of 
space where the complex field exhibits high peak power 
fluctuations at the initial time t = 0, see e. g. the region 
where x G [-50,-40] in Fig. 2 and in Fig. 3(a),(c),(e). 
While the localized peaks shown in Fig. 2(a) drift with 
small velocities in the (x, t) plane, their peak power in¬ 
creases with time. This increase of the peak power of 
the localized structures goes simultaneously with a re¬ 
duction of their spatial width. The typical spatial scale 
of the power fluctuations evolves from a value of ^ 10 at 
initial stage {t = 0) to the healing length of ^ 1 at long 
evolution time {t > 10) (see Fig. 3(a) , (c), (e)). 

A clear signature of the change in the fluctuation scale 
can be observed in the Fourier space. As shown in Fig. 
3(b), 3(d) and 3(f), the power spectrum |'^(fc,t)p of the 
random wave is indeed found to significantly broaden 
with time. Let us emphasize that the power spectrum 











5 


of the wave broadens while always keeping exponentially 
decaying wings. 



Figure 2. Numerical simulations of Eq. m in focusing 
regime {a = +1). (a) Spatio-temporal evolution of the power 
lipix, t)|^ of the wave while starting from the random complex 
field having a gaussian power spectrum and a gaussian statis¬ 
tics, see Fig. 1. (b) Corresponding time evolution of linear 
(kinetic) energy Hl and of nonlinear energy Hnl- 


The evolution of the typical space scales and peak pow¬ 
ers of the random fluctuations that is described above 
goes together with a process of energy balance between 
linear and nonlinear effects. Fig. 2(b) shows the evo¬ 
lutions in time of the linear (kinetic) energy and of 
the nonlinear energy that are associated with the 

spatio-temporal evolution plotted in Fig. 2(a). The wave 
system is initially placed in a highly nonlinear regime in 
which the nonlinear energy is one order of magnitude 
greater than the linear energy (|i?ArL| ~ 10\Hl\). As 
a result of nonlinear propagation, linear and nonlinear 
effects come into balance and after a short transient evo¬ 
lution, the wave system reaches a state in which linear 
and nonlinear energies have the same order of magnitude 
(I iJvTL I 2 1 Hi I, see Fig. 2(b) ). 

It is interesting to compare the spatio-temporal evo¬ 
lution shown in Fig. 2(a) with the one shown in Fig. 


2(a) of ref. 


In ref. [S^, the authors study the 


nonlinear evolution in space and time of a condensate 
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Figure 3. Numerical simulations of Eq. m in focusing regime 
(cr = -1-1). Spatial evolution of the power |'i/i(x, t)f of the ran¬ 
dom wave at times (a) ti = 0.24, (c) ti — 2.44 , (e) = 16.04. 

Power spectra (red lines) of the random wave at 

times (b) fi = 0.24, (d) ti = 2.44 , (f) fi = 16.04. The spectra 
plotted in blue lines represent the gaussian power spectrum 
of the random initial complex field defined by Eq. m and 
Eq. 0. 


perturbated by a small noise, i. e. ipoi^) = 1 + 
where |??(a;)| << 1 is a small noise with broad spectrum. 


As shown in Fig. 2(a) of ref. [82|, the fact that there is 


only a weak random modulation of the initial condition 
gives rise to a spatiotemporal diagram in which localized 
structures are distributed in space and time in a way 
that is more regular than the random pattern shown in 
Fig. 2(a) of this paper. Starting from our initial con¬ 
dition with a broad gaussian spectrum, there are some 
wide regions of space in which no localized structures are 
observable while there are some other regions of space in¬ 
cluding many localized structures. As shown in ref. [s^, 
the localized structures emerging from a complex field 
having initially a gaussian statistics can be locally fitted 
by some analytical functions corresponding to solitons 
on finite background, such as e. g. the Peregrine soliton. 
Fitting procedures implemented in the numerical work 
presented in ref. have shown that many solitons on 
finite background can be also found while seeding the 
wave system from a condensate perturbated by a small 
noise. 


Despite localized structures looking like solitons on 
finite background can be observed while starting from 
those two different random initial conditions, signifi¬ 
cantly different statistical features are observed at long 
evolution time. It has been shown in ref. Q that gaus¬ 
sian statistics emerges from the nonlinear evolution of the 
noisy condensate. As shown in Fig.4(a), heavy-tailed de¬ 
viations from gaussianity are contrarily found to emerge 
from a random complex field having initially a gaussian 
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statistics [s^ . It is an open question to understand how 
the interplay among localized structures does not pro¬ 
duce the same statistics at long evolution time while 
starting from different noisy initial conditions. 


Rational solutions of the ID-NLSE such as Akhme- 
diev breathers are considered as prototype of rogue waves 
Ulil. These coherent structures have been generated 
in optical fiber experiments [^, 761 and in hydrodynami- 
cal experiments |72l ] . The interaction of solitons |19l| and 
the collision of breathers haye been studied theoretically 
and experimentally 7^. Our work points out the rel- 
eyance of these works and the need to extend it in the 
context of random nonlinear wayes. In particular, it is 
an open question to understand the mechanisms of the 
emergence of coherent structures from the two different 
initial conditions, i.e. the plane wave with small noise on 
one hand and the random wave computed from the RP 
model on the other hand. 


Fig. 4 shows the time evolution of the PDF of 
|'0(a;,t)p. The nonlinear random field has a statistical 
evolution in which the PDF of power fluctuations con¬ 
tinuously moves from the exponential distribution (plot¬ 
ted in red line in Fig. 4) to the heavy-tailed distribution 
plotted in magenta line in Fig. 4. For t > 10, the wave 
system reaches a statistical stationary state in which the 
PDF no longer changes with time [SOj. This statistical 
stationary state is determined by the interaction of coher¬ 
ent nonlinear structures such as for instance Akhmediey 
breathers, Kuznetsoy-Ma solitons and also linear disper- 
siye radiation [IS [oS [sHI . It is now an open question to 
determine the mechanisms in integrable turbulence that 
lead to the establishment of a stationary state in which 
statistical properties of the waye system do not change 
in time. Tools from the inyerse scattering theory could 
be used to investigate this question of fundamental im¬ 
portance 17, 95-97, 9^ 



Figure 4. Numerical simulations of Eq. m in focusing regime 
(a = -1-1). PDF of the power at time f = 0 (red line), 

t = 0.2 (green line), t = 0.4 (blue line), f = 10 and t* = 20 
(magenta line). The PDF is stationary from t ~ 10, i.e. the 
PDF plotted in magenta line does not change with time for 
t > 10. 


Defocusing regime 


In this Section, we consider the defocusing regime 
(ct = —1) and we use numerical simulations of Fq. ([T]) to 
investigate the propagation of a partially coherent wave 
generated at t = 0 by a random complex field identical 
to the one used in Sec. □ 



Figure 5. Numerical simulations of Eq. o in defocusing 
regime (a = —1). (a) Spatio-temporal evolution of the power 
\'ip{x,t)f of the wave while starting from the random complex 
field having a gaussian power spectrum and a gaussian statis¬ 
tics, see Fig. 1. (b) Corresponding time evolution of linear 
(kinetic) energy Hl and of nonlinear energy Hnl- 


Fig. 5(a) shows the spatio-temporal evolution of the 
partially-coherent wave seeded by a random initial field 
identical to the one used in Fig. 2(a). Spatiotemporal 
features emerging from the nonlinear propagation in de¬ 
focusing regime drastically contrast with those found in 
the focusing regime. Instead of bright localized struc¬ 
tures, we now observe the emergence of dark localized 
structures propagating at various speeds in the {x, t) 
plane. Fig. 6(a) and 6(b) show that the initial stage 
of nonlinear evolution is now characterized by a fast de¬ 
cay of the peaks of highest intensities, see e. g. the region 
where x S [—50, —40]. During the initial evolution of the 
random wave, the leading and trailing edges of peaks of 
highest intensities strongly sharpen. This leads to some 
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gradient catastrophes which are regularized by the gen¬ 
eration of dispersive shock waves (DSWs) [ 73 , lUii. 
As in the focusing case, the typical spatial scale of the 
random fluctuations decreases from ~ 10 at the initial 
stage {t = 0) to the healing length of ^ 1 at long evo¬ 
lution time {t > 10). The stochastic evolution shown in 
Fig. 6(c) is determined by the interaction of nonlinear 
coherent structures such as dark solitons or DSWs and 
of linear radiation. 

Fig. 6(b), 6(d) and 6(f) show that nonlinear propaga¬ 
tion in defocusing regime induces a spectral broadening of 
the random wave. This spectral broadening phenomenon 
is quantitatively less pronounced than the one observed 
in the focusing regime, see Fig. 3(b), 3(d), 3(f). However 
the power spectrum of the wave broadens while always 
keeping exponentially decaying wings, as in the focusing 
regime. 
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Figure 6. Numerical simulations of Eq. o in defocusing 
regime (a = —1). Spatial evolution of the power \-ijj{x,t)\'^ of 
the random wave at times (a) ti = 0.24, (c) = 2.44 , (e) 

ti — 16.04. Power spectra (red lines) \'ip{k,t)\‘^ of the random 
wave at times (b) ti — 0.24, (d) ti = 2.44 , (f) ti — 16.04. 
The spectra plotted in blue lines represent the gaussian power 
spectrum of the random initial complex field dehned by Eq. 
@ and Eq. (|4]). 


5(b). 



Figure 7. Numerical simulations of Eq. 0 in defocusing 
regime (cr = —1). PDF of the power at time t — 0 

(red line), t — 0.2 (green line), t = 1 (blue line), f = 10 and 
t* = 20 (magenta line). The PDF is stationary from t ~ 10, 
i.e. the PDF plotted in magenta line does not change with 
time for t > 10. 

In the defocusing regime, the nonlinear random field 
has a statistical evolution in which the PDF of power 
fluctuations continuously moves from the exponential dis¬ 
tribution (plotted in red line in Fig. 7) to the low-tailed 
distribution plotted in magenta line in Fig. 7. As in 
the focusing regime, the wave system exhibits a statis¬ 
tical stationary state and the PDF computed at t = 10 
(magenta line in Fig. 7) does not change anymore with 
time 0. This statistical stationary state is determined 
by the interaction of coherent nonlinear structures such 
as for instance dark solitons, dispersive shock waves and 
also linear radiation. As for the focusing regime, it is an 
open question to determine the mechanisms in integrable 
turbulence that lead to the establishment of a stationary 
state in which statistical properties of the wave system 
do not change in time. Tools of the inverse scattering 
transform could be of interest for the investigation of this 
question [^ . 


OPTICAL FIBER EXPERIMENTS IN FOCUSING 
AND IN DEFOCUSING PROPAGATION 
REGIMES 


The spatioteinporal evolution shown in Fig. 5(a) and 
in Fig. 6(a), (c), (e) goes together with a process of en¬ 
ergy balance between linear and nonlinear effects. Fig. 
5(b) shows the time evolutions of the linear (kinetic) en¬ 
ergy Hl and of the nonlinear energy Hi\fL that are as¬ 
sociated with the spatio-temporal evolution plotted in 
Fig. 5(a). The wave system is initially placed in a 
highly nonlinear regime in which the nonlinear energy 
is one order of magnitude greater than the linear energy 
{Hnl ^ IOFTl). As in the focusing regime, linear and 
nonlinear effects come into balance and after a short tran¬ 
sient evolution, the wave system reaches a state in which 
Hn]^ and have the same order of magnitude, see Fig. 


Optical hber experiments provide versatile and pow¬ 
erful tabletop laboratory to investigate the complex dy¬ 
namics of IDNLSE, hydrodynamic-like phenomena and 
the statistical prop erties of n onlinear random waves 

[1 i, US El S EJ 0, H Eiil. 


One of the most critical constraint of these experiments 
is the finite spectral bandwidth of usual detectors. The 
typical response time of the fastest detector and oscillo¬ 
scope is several tens of picoseconds. On the other hand, 
with the usual parameters of standard experiments using 
optical fibers, the typical “healing time” scale character¬ 
izing the equilibrium between the nonlinearity and the 
dispersion is around one picosecond [s^. For this rea- 





































son, the picosecond is also the order of magnitude of the 
time scale associated to the power fluctuations of par¬ 
tially coherent fiber lasers [H y, [13 ■ 

As a consequence spectral filters are therefore often 
used to reveal extreme events in time-domain experi¬ 
ments d, 23, 10, In the case of pulsed experiments, 
it is possible to evidence shot-to-shot spectrum fluctua¬ 
tions with a dispersive Fourier transform measurement 
m [ 23 , [n m. To the best of our knowledge, up to 
our recent works [0,1^ ; accurate and well-calibrated 
measurement of the PDF characterizing temporal fluc¬ 
tuations of the power of random light with time scale of 
the order of picosecond had never been performed. 
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X = 457nm 


Sampled signal 


signal at a wavelength A = 457nm. 

In our experiments, the linearly polarized partially co¬ 
herent wave is emitted by a “continuous” wave (cw) Yt¬ 
terbium fiber laser at As = 1064nm. This cw laser 
emits numerous (typically 10 ^) uncorrelated longitudi¬ 
nal modes. The reader can refer to [s^ for the details of 
the experimental setup and of the statistics measurement 
procedure. 

In this paper we present results obtained with two 
different fibers (hber 1 and 2 ) having opposite sign of 
the group velocity dispersion at the wavelength of the 
signal As = 1064nm. The fiber 1 is a 15m-long highly 
nonlinear photonic crystal fiber (provided by Draka 
France company) with a nonlinear third order coefficient 
7 ~ 50W“^km“^ and a group velocity dispersion 

coefficient (32 — —2Qps^/km. The fiber 2 is a lOOm-long 
polarization maintaining fiber with a nonlinear third 
order coefhcient 7 ~ 6 W“^km“^ and a group velocity 
dispersion coefficient /?2 — -\-2t)ps^/km. We launch 

a mean power < P >= 0.6W in the experiments 
performed with the fiber 1 (focusing case) and a mean 
power < P >= 4.W in the experiments performed 
with the fiber 2 (defocusing case). Note that the results 
obtained with the fiber 1 have been presented in detail in 
whereas the results obtained with the fiber 2 are new. 


b 



Figure 8. Experimental measurement of the statis¬ 
tics of random light a. Principle. Optical sampling of 
the partially-coherent wave fluctuating with time (the signal) 
is achieved from sum frequency generation in a second order 
(x^^^) crystal. Blue pulses are generated at A = 457nm from 
the interaction of the signal with periodic femtosecond pump 
pulses inside a crystal. PDF is computed from the peak 
powers of the blue pulses, b. Nonlinear propagation in 
optical fiber The initial partially coherent light is launched 
inside a single mode optical fiber (either in the focusing (fiber 
1) or defocusing (fiber 2) regime of dispersion. The statistics 
and the spectrum is measured before and after the propaga¬ 
tion in the fiber. 

We have developed an original setup based on asyn¬ 
chronous optical sampling (OS) which allows the precise 
measurement of statistics of random light rapidly fluctu¬ 
ating with time (see Fig. [El). 

The “random light” under investigation is a partially 
coherent wave and it is called the “signal”. The signal is 
optically sampled with 140fs-pulses and the PDF is com¬ 
puted from the samples. The optical sampling is obtained 
from the second order nonlinearity in n BBO crys¬ 
tal. The sum-frequency generation (SFG) between the 
signal at Ag = 1064nm and short “pump” pulses having 
a central wavelength Ap = 800nm provide samples of the 


We first measure the PDF at the output of the laser. 
In all experiments presented in this letter, the mean out¬ 
put power of the Ytterbium laser is fixed at (P) = lOW. 
At this operating point, the statistics of the partially co¬ 
herent wave follows the normal law. Indeed, as plotted in 
red in Figs. Ill c and Eld, the PDF of the normalized power 
P/ (P) is very close to the exponential function. Assum¬ 
ing that the real part and the imaginary parts are statisti¬ 
cally independent, this exponential distribution of power 
corresponds to a gaussian statistics of the field. It is im¬ 
portant to note that the dashed black lines in Figs. Elc 
and Eld are not a fitted exponential function but repre¬ 
sent the exact normalized PDF[P/{P)] = exp(—P/(P)). 
To the best of our knowledge, PDF of so rapidly fluctu¬ 
ating optical signals has never been quantitatively com¬ 
pared to the normalized exponential distribution. 

We use the output of the laser as a random source and 
we launch the partially coherent signal into optical fibers 
1 and 2. Experiments have been carefully designed to 
be very well described by the ID-NLSE. In particular, 
the signal wavelength As = 1064nm is far from the 
zero-dispersion wavelength (Ag— 970nm for the fiber 1 
and Ao > 1300nm for the hber 2). Moreover the optical 
spectral widths (see Figs. Ela and Elb) remain suffi¬ 
ciently narrow to neglect stimulated Raman scattering 
(SRS) and high-order dispersion effects. The linear 
losses experienced by optical fields in single pass in 
the hbers are neglictible. These total losses are around 
0.3% in the fiber 1 and around around 2.5% in the fiber 2. 
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Figure 9. Experiments. a,b : Optical spectra. In¬ 
put spectrum (red line). Spectrum at the output of fiber 
1 (green line) and fiber 2 (blue line). The circles represents 
optical spectra computed from numerical simulation of ID- 
NLSE. c, d : PDF of normalized optical power Pj(P) plot¬ 
ted in logarithmic scale. Normalized exponential distribution 
PDF[P/ < P >] = exp(—P/ < P >) (black dashed line). 
PDF of the input random light (red line). PDF of the light at 
the output of the fiber 1 (green solid line) and at the output 
of the fiber 2 (blue solid line). The dashed lines represent the 
PDF computed from numerical simulation of ID-NLSE with 
the parameters of the experiments. The inset represents the 
same PDF in linear scale. 


Despite the broadening of the optical spectrum is of 
nearly the same importance in focusing and in defocus- 
ing regimes (see Figs. Oa andOb) our experiments reveal 
that the distortions of the statistics of the random waves 
strongly depend on the sign of the group velocity disper¬ 
sion coefficient. 

The experiments performed in the focusing regime 
(fiber 1) reveal the occurrence of numerous extreme 
events (RW) (see green curve in Fig. Oc). The com¬ 
parison between the initial PDF (see red line in Fig. Oc) 
and the output PDF (see green curve Oc) shows an im¬ 
pressive change in the statistical distribution of optical 
power. The initial field follows the normal law and its 
PDF is an exponential function whereas the output PDF 
of optical power exhibits a strong heavy-tail. 

On the other hand, the PDF experimentally measured 
at the output of fiber 2 in the defocusing regime exhibits 
a very low tail (see Fig. |9ld). Light fluctuations of a 
high power are found with a probability that has been 
strongly reduced as compared to the normal law. More¬ 
over, contrary to the initial exponential distribution, the 
most probable value for the power is not the zero value 
(see inset of Fig. Id) ). 

Numerical simulations show that experiments pre¬ 
sented above are very well described by the integrable 
ID-NLSE. The initial conditions are computed from the 
random phase assumption as in the section „ We have 


performed Monte Carlo simulations with ensemble av¬ 
erage over thousands of realizations. We integrate the 
ID-NLSE with experimental parameters : 



^2 

YW 


7|i/;|V 


(5) 


where (32 is the group velocity dispersion coefficient 
and 7 is the effective Kerr coefficient. Optical spectra 
and PDEs computed from the numerical integration of 
the ID-NLSE are in quantitative agreement with exper¬ 
iments both in the focusing and in the defocusing cases 
(see dashed green and blue lines in Fig. IH]). 

Moreover, the numerical simulations show that inte¬ 
grable turbulence is characterized by a statistical station¬ 
ary state both in focusing and defocusing regime (see Fig. 
[ini)- In particular, the average of the nonlinear and lin¬ 
ear parts of the Hamiltonian ((Hatl) and {Hl}) evolves 
to constant values (see Figs. [TOla and ITOlb). Note that 
we represent here the average of Hnl and Hl over hun¬ 
dreds of realizations whereas the values oi Hnl and Hl 
computed on only one realization are plotted in section □ 

Fig. llOl c represents PDFs computed from numerical 
simulations for different lengths of propagation in the fo¬ 
cusing regime of dispersion. The red line is the PDF 
at z = Om, the green line is the PDF at z = 15m cor¬ 
responding to the experiments (see Fig. Elc) and the 
black line corresponds at the stationary PDF (computed 
at z = 500m). Fig. [TUld represents PDFs computed from 
numerical simulations for different lengths of propagation 
in the defocusing regime of dispersion. The red line is the 
PDF at z = Om, the blue line is the PDF at z = 100m 
corresponding to the experiments (see Fig. Od) and the 
black line corresponds to the stationary PDF (computed 
at z = 500m). 

Note that comparable deviations from gaussian statis¬ 
tics have been reported in ID “spatial experiments” in 
which the transverse intensi ty p rofile of optical beams 
randomly fluctuates in space j22l| . In these spatial exper¬ 
iments performed in focusing and defocusing regime, the 
speckle fields are localized and random waves decay to 
zero at infinity [iSlii. 1ST with usual zero-boundary 
conditions has been used in ref. 0 to describe these ex¬ 
periments. In the long-term evolution of the wave system 
with zero boundary conditions, solitons separate from 
dispersive waves in the focusing regime. In the defocus¬ 
ing regime, only dispersive waves persist at long evolution 
time. 

On the contrary, our experiments and numerical simu¬ 
lations are made with non-zero boundary conditions and 
with non-localized random waves. This widens the per¬ 
spectives of experimental integrable turbulence studies. 
With random waves having an infinite spatial extension, 
solitons and dispersive waves never separate from each 
other and they always interact. Moreover breathers and 
solitons on finite background can emerge from nonlinear 
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a b 






Figure 10. Numerical Simulations, a and b : evolution 
of the average of the nonlinear Hnl and linear Hl Hamilto¬ 
nians in the focusing (a) and defocusing (b) cases. Hnl and 
Hl are normalized to the value of the Hamiltonian (constant 
of motion) H = Hnl + Hl before the averaging over hun¬ 
dreds of realizations, c. PDFs computed at different 2 
in the focusing case. 2 = Om (red line), 2 = 15m (green 
line, identical to dashed green line in Fig. |^c) and 2 = 100m 
(stationary PDF), d. PDFs computed at different 2 in 
the defocusing case. 2 = Om (red line), 2 = 100m (blue 
line, identical to dashed blue line in Fig. [9]d) and 2 = 500m 
(stationary PDF). 


interaction in the focusing regime. In the defocusing case, 
the fact that the random field does not decay at infinity 
means that dark solitons can be sustained and interact 
all-together with dispersive waves at any time (any value 
of 2 in onr fiber experiments). 

As a conclusion of Sec. „ we have experimentally stud¬ 
ied the evolution of the statistics of random waves whose 
propagation is very well described by ID-NLSE both in 
focusing and defocusing case. 


In the defocusing case, we have experimentally and 
numerically demonstrated that the probability of occur¬ 
rence of large waves decreases as a result of nonlinear 
propagation. 


In the focusing case, we have evidenced the statistical 
emergence of RW from nonlinear propagation of random 
light. In [^, we have also shown that solitons on fi¬ 
nite background such as Akhmediev breathers, Peregrine 
solitons or Kuznetsov-Ma solitons having a short dura¬ 
tion and a high power seem to emerge on the top of the 
highest fluctuations. This strengthens the idea that the 
emergence of deterministic solutions of ID-NLSE such 
as Akhmediev breathers in nonlinear random fields is 
a major mechanism for the formation of rogue waves 
(5I. [itt 6^. Note that the emergence of such co¬ 

herent structures in incoherent fields has been already 
theoretically studied in non integrable wave turbulence 
0 [zi and in integrable turbulence emerging from a 
modulationaly unstable condensate [^, |4l|, [0 . 


Note finally that in one dimensional deep water exper¬ 
iments, relatively small deviations from gaussianity have 
been observed and interpreted in the framework of wave 
turbulence theory H, 5^. On the contrary, our op¬ 
tical fiber setup provides an accurate laboratory for the 
exploration of strongly nonlinear random wave systems 
ruled by the ID-NLSE. 


SEPARATION OF SCALES AND 
INTERMITTENCY PHENOMENON 

Statistical features presented in Sec. □ and in Sec. □ are 
relevant to global random fields in the sense that all the 
fluctuations scales of the random waves are taken into 
account and contribute to the statistics. However, sepa¬ 
rating large scales from small scales is known to provide 
rich statistical information about nonlinear systems of 
random waves. In this respect, the phenomenon of inter- 
mittency is defined in the general context of turbulence 
as a departure from the Gaussian statistics that grows 
increasingly from large scales to small scales [s^. 

Following the definition given by Frisch in ref. [83| . 
a random function R(x) of space x is defined as being 
intermittent when it displays some activity over a frac¬ 
tion of space that decreases with the scale under consid¬ 
eration. Considering stationary random processes, the 
intermittency phenomenon can be evidenced and quanti¬ 
fied by using spectral filtering methods. The existence of 
deviations from gaussianity is usually made through the 
measurement of the kurtosis of the fluctuations that are 
found at the output of some frequency filter. Considering 
the high-pass filtered signal (x) defined in the spatial 
domain as 

R{x) = J R{k), (6) 

R>{x)= [ dke^^^R{k), (7) 

J\k\>i 

the random function R(x) is intermittent at small scales 
if the kurtosis: 


<e) = 


m{x)w 


( 8 ) 


grows without bound with the filter frequency ^ [83| . 

Although we are going to use here the definition of in¬ 
termittency given by Frisch (Eqs. (I6])-(I8|)), it should be 
emphasized that the exact nature of the spectral hltering 
process is not of a fundamental importance and that the 
intermittency phenomenon can also be evidenced from 
the use of various frequency Liters. Spectral fluctuations 
can be examined at the output of an ideal one-mode 
spectral filter passing only a single Fourier component 
0, 0. Time fluctuations at the output of bandpass 
frequency Liters can also be considered [sj, l83|,0. PDFs 
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of second-order differences of the wave height have also 
been measured in wave turbulence [85| . Using this kind 
of filtering techniques, the phenomenon of intermittency 
has been initially reported in fully developed turbulence 

S I but it is also known to occur in wave turbulence jS^ 
, solar wind (89| or in the Faraday experiment [87l| . 
So far, the intermittency phenomenon has been ascribed 
to physical systems that are described by non-integrable 
equations. We are going to use numerical simulations 
of Eq. ([T]) to show that intermittency is a statistical 
phenomenon that also occurs in the field of integrable 
turbulence Q. 



Figure 11. Numerical simulations of Eq. m in the focus¬ 
ing regime {a = -|-1) for t* = 20. (a) Spatial fluctuations 
of R{x)/ < R{x) >= 5R('0(a;, t*))/ < di{ilj{x,t*)) > that are 
found at the output of an ideal highpass frequency filter hav¬ 
ing a cutoff frequency ^ = 0 (Eq. ([6])) and (d) corresponding 
PDF of R{x)/ < R{x) >. (b) and (e), same as (a) and (d) 
for ^ = 1. (c) and (f), same as (a) and (d) for ^ = 2. The 
insets in (a), (b), (c) represents the Fourier power spectra of 
the random fields plotted in (a), (b), (c). 


Taking the wave system extensively described in Sec. „ 
we now consider the spatial fluctuations of the real part 
t*)) of the complex field ^(x, t*) that has reached 
the stationary statistical state at t* = 20. In other words, 
the filtering process and the statistical treatment defined 
by Eqs. ©-(0) are now applied to the random vari¬ 
able R{x) = 5R('0(x, t*)). Spatial and statistical features 
found at the output of the highpass frequency filter are 
shown in Fig. 11 for the focusing regime. 

When ^ = 0, the random process R{x) is not filtered 
and Fig. 11(a) shows the spatial evolution of the random 
variable R{x)/ < R{x) >= ^{'ij}{x,t*))/ < 5i(^(x,f*)) > 
in this situation. Fig. 11(d) represents the correspond¬ 
ing PDF of R{x)/ < R{x) >. Without any filtering pro¬ 
cess, the heavy-tailed deviations from gaussianity shown 


in Fig. 11(d) are identical to those already shown in Fig. 
4 but for |'0(x,t*)p. Increasing the cutoff frequency ^ of 
the ideal highpass filter (Eq.©), fluctuations of smaller 
and smaller scales are observed at the output of the high- 
pass filter together with peaks of higher and higher am¬ 
plitudes, see Fig. 11(b) and 11(c). Fig. 11(e) and Fig. 
11 (f) show that deviations from gaussianity become heav¬ 
ier when the cutoff frequency ^ of the highpass filter is 
increased. These statistical features represent qualitative 
signatures of the intermittency phenomenon. Comput¬ 
ing the kurtosis k(^) of R{x) = 5J('0(x, t*)) for increasing 
values of the cutoff frequency we find a monotonic 
increase that complies with the definition of the inter¬ 
mittency phenomenon given by Frisch, see Fig. 13(a) 

M- 



Figure 12. Numerical simulations of Eq. o in the defocus- 
ing regime {a = —1) for t* = 20. (a) Spatial fluctuations 
of R(x)/ < R(x) >= ?R.(ip{x,t*))/ < K('i/'(x, t*)) > that are 
found at the output of an ideal highpass frequency filter hav¬ 
ing a cutoff frequency ^ = 0 (Eq. ©) and (d) corresponding 
PDF of R{x)/ < R{x) >. (b) and (e), same as (a) and (d) 
for ^ = 1. (c) and (f), same as (a) and (d) for ^ = 2. The 
insets in (a), (b), (c) represents the Fourier power spectra of 
the random fields plotted in (a), (b), (c). 


As shown in Fig. 12, features qualitatively similar 
to those described for the focusing regime are found 
in the defocusing regime. Fig. 12(a) (resp. Fig. 12 
(d)) shows the spatial evolution (resp. the PDF) of 
R{x)/ < R{x) >= 5i('0(x,t*))/ < K('0(x,t*)) > for 
^ = 0, when the random process is not hltered. The low¬ 
tailed deviations from gaussianity shown in Fig. 12(d) 
are identical to those already shown in Fig. 7 but for 
|'0(x,t*)p. They characterize the stationary statistical 
state found in the defocusing regime. Increasing the cut¬ 
off frequency ^ of the ideal highpass filter (Eq.©), fluc¬ 
tuations of smaller and smaller scales are observed at 
the output of the highpass filter together with peaks of 
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higher and higher amplitudes, see Fig. 12(b) and 12(c). 
Fig. 12(e) and Fig. 12(f) show that deviations from PDF 
computed for ^ = 0 become heavier when the cutoff fre¬ 
quency ^ of the highpass filter is increased. As shown 
in Fig. 13(b), the kurtosis k(^) monotonically increases 
with as for the focusing regime. Let us recall that a 
kurtosis k equal to 3 corresponds to a random field having 
a gaussian statistics. The fact that the initial value of the 
kurtosis k(^ = 0) is lower (resp. greater) than 3 in defo- 
cusing (resp. focusing) regime complies with the fact that 
the unfiltered field exhibit low-tailed (resp. heavy-tailed) 
deviations from gaussianity at t* = 20. Note that the 
growth of the kurtosis k{^) is greater in focusing regime 
than in defocusing regime, see Fig. 13. 



Figure 13. Numerical simulations of Eq. (IT|) for t* = 20. Evo¬ 
lution of the kurtosis k as a function of the cutoff frequency ^ 
of the highpass filter defined by Eq. ®. (a) Focusing regime 
{a = -1-1). (b) Defocusing regime (cr = —!)• 


CONCLUSION 

The work presented in this paper deals with the general 
question of statistical changes experienced by ensembles 
of nonlinear random waves propagating in systems ruled 
by integrable equations. It enters within the framework 
of “integrable turbulence” which is a new field of research 
introduced by Zakharov to address specifically this ques¬ 
tion that “composes a new chapter of the theory of tur¬ 
bulence” [i^. In our work, we have specifically focused 
on optical fiber systems accurately described by the inte¬ 
grable one-dimensional nonlinear Schrodinger equation. 
We have considered random complex fields having a gaus¬ 
sian statistics and an infinite extension at initial stage. 
Numerical simulations with periodic boundary conditions 
and optical fiber experiments have been used to investi¬ 
gate spectral and statistical changes experienced by non¬ 


linear waves both in focusing and in defocusing propaga¬ 
tion regimes. 

As a result of propagation in the strongly nonlinear 
regime, the power spectrum of the random waves is found 
to broaden while taking exponential wings both in focus¬ 
ing and in defocusing regimes, see Sec. □ and Sec. □ 
In the nonlinear regime, this spectral broadening phe¬ 
nomenon is a signature of a process in which the typical 
spatial scale of the random fluctuations decreases to reach 
the healing length of the wave system at long evolution 
time. Numerical simulations have revealed that heavy¬ 
tailed deviations from gaussian statistics occur in the fo¬ 
cusing regime while low-tailed deviations from gaussian 
statistics are found in the defocusing regime. These sta¬ 
tistical behaviors have been observed in optical fiber ex¬ 
periments relying on the implementation of an original 
and fast detection scheme, see Sec. „ Numerical sim¬ 
ulations made at long evolution times have also shown 
that the wave system exhibits a statistically stationary 
state in which neither the PDF of the wave field nor the 
spectrum change with the evolution variable. Separat¬ 
ing fluctuations of small scale from fluctuations of large 
scale, we have finally revealed the phenomenon of in- 
termittency; i.e., small scales are characterized by large 
heavy-tailed deviations from Gaussian statistics, while 
the large ones are almost Gaussian. This intermittency 
phenomenon has been observed both in focusing and in 
defocusing propagation regimes, see Sec. □ 

As underlined in Sec. „ the determination of the PDF 
of a random wave field represents an issue of importance 
in the field of integrable turbulence. Given a nonlin¬ 
ear partial differential integrable equation together with 
some given initial and boundary conditions, there is no 
systematic theory that allows one to determine the PDF 
of the wave field at asymptotic stage (i.e. long evolution 
time). 

The wave turbulence (WT) theory describes out-of¬ 
equilibrium statistical mechanics of random nonlinear 
waves in the weakly nonlinear regime The WT 

theory is commonly used to treat wave systems that 
are ruled by non integrable equations and in which the 
long-time evolution is dominated by resonant interactions 
among waves {i.e. Fourier components). Using the clo¬ 
sure of the statistical moments of the random field and 
neglecting non resonant interactions, kinetic equations 
describing the long-term evolution of the wave spectrum 
can be derived. However it has been shown that the short 
term e volut ion of WT can be influenced by non-resonant 
terms 124l |. In the specific case of the integrable ID- 
NLSE, there are only trivial resonance conditions 
As a result, collision terms found in the kinetic equations 
determined from the WT theory vanish [13, ■ Us¬ 

ing the closure of the moments and keeping the contri¬ 
bution of non resonant terms, it is however possible to 
derive “quasi-kinetic” equations th at de scribe the evolu¬ 
tion of the wave spectrum 14, 42, I23l| . This approach 
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has been used to study the evolution of the statistics both 
in focusing and defocusing case in ref. 4^. This kind 
of WT treatment has been shown to properly describe 
the evolution of the kurtosis of the wave field together 
with spectral c hang es that occur in the weakly nonlinear 
regime 

There are now many open questions regarding inte- 
grable turbulence in the strongly nonlinear regime dis¬ 
cussed throughout this paper. The inverse scattering the¬ 
ory (1ST) provides a natural framework for the investi¬ 
gation of statistical properties of nonlinear wave systems 
ruled by integrable equations. In particular, 1ST has been 
used in ref. 17| to describe some experiments examining 


nonlinear diffraction of localized incoherent light beams 
[ 2 ^ . The theoretical description of this experiment can 
be made by using standard tools of the 1ST because it is 
fully compatible with the central assumption of 1ST that 
the wave field decays at infinity. 1ST has thus been used 
both in the focusing and in the defocusing regimes to de¬ 
termine some mathematical expressions for the PDF of 
the wave field [ 3 . 

Our experiments and numerical simulations made with 
non-decaying and non-localized random waves open new 
questions about statistical properties of incoherent waves 
in integrable turbulence. First, the statistical properties 
characterizing the asymptotic stage of integrable turbu¬ 
lence are of fundamentally different natures for waves 
fluctuating around a constant background and for decay¬ 
ing waves. With continuous random waves having an 
infinite spatial extension, the nonlinear evolution of the 
random wave can no longer produce individualized soli- 
tons at long time. On the other hand, solitons never sep¬ 
arate from each other and they always interact. Moreover 
solitons on finite background can emerge from nonlinear 
interaction in the focusing regime. In the defocusing case, 
the fact that the random field does not decay at infinity 
means that dark solitons can be sustained and interact 
at any time (any value of z in our fiber experiments). 

In our work, the relevant boundary conditions are pe¬ 
riodic boundary conditions in a box of size L. Random 
waves of infinite spatial extension such as the ones con¬ 
sidered in experiments reported in Sec. □ can be described 
by taking the limit of a box of infinite spatial extension 
{L —>• 00 ). Rigorously speaking, the theoretical frame¬ 
work for dealing with integrable wave systems and peri¬ 
odic boundary conditions is finite gap theory |llOl |. So 
far, no theoretical work has been made in this framework 
to determine statistical properties of ensembles of nonlin¬ 
ear random waves. However some recent results point out 
the possibility to use the inverse scattering transform for 
the focusing nonlinear Schrbdinger equation with nonzero 
boundary conditions |l 22 | . 

The gaussian statistics of the initial condition is a key 
point of our experimental work. In the focusing regime 
of the ID-NLSE, the statistics of the field measured in 
the statistically stationary state strongly depends on the 


nature of the initial condition. In our experiments and 
numerical simulations made with a complex field hav¬ 
ing initially a gaussian statistics, heavy-tailed deviations 
from gaussianity have been observed in the statistically 
stationary state (i.e. at long time), as discussed in Sec. 
n and in Sec. □ If the initial condition is now made from 
a plane wave (or a condensate) with an additional noise, 
the nonlinear stage of modulational instability is charac¬ 
terized by a stationary statistics following the normal law 
0. The fact that there exists such a strong qualitative 
difference between statistics measured in the stationary 
state while starting from different noisy initial condition 
is an intriguing issue. 

We hope that our results will stimulate new theoretical 
works aiming to understand the mechanisms leading to 
the strongly non gaussian statistics in focusing and defo¬ 
cusing ID-NLSE systems with random initial conditions 
and non-zero boundary conditions. 
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